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The origins of recently reported anisotropy of the local velocity field of nearby galaxies (velocities < 500 
1 km/s corresponding to the distance less than 8 Mpc) are studied. The exact solution of the Newtonian 



equation for the expanding Universe is obtained. This solution allows us to separate the Newtonian motion 
of nearby galaxies from the Hubble flow by the transition to the conformal coordinates. The relation between 
the Hubble flow and the Newtonian motion is established. We show that the anisotropic local velocity field 
of nearby galaxies can be formed by such a Newtonian motion in the expanding Universe, if at the moment 
of the capture of galaxies by the central gravitational field. 
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I. INTRODUCTION 



Recent observation of the local velocity field of galaxies gives a three-dimensional ellipsoid with 
different values of the Hubble parameter, clearly showing its anisotropic character 0, |^ . 

In this paper, we present a possible point of view that this local velocity field of galaxies can be 
explained by their Newtonian motions. 

The analysis of the observational data will be based on the radial velocities of nearby galaxies, 
belonging to the Local Group. Our paper is organized in the following manner. In Section II, the 
cosmic evolution is described. In Section III, we introduce the Newtonian motion and separate this 
motion from the cosmic one. In Section IV, the initial data of the galaxies capture by a central 
gravitational field is considered. In Sections V, the simplest example is given to elucidate our results. 
The paper ends with the conclusions. 



II. COSMIC EVOLUTION OF A FREE PARTICLE 



Effects of cosmic evolutions are considered in the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metrics 

(ds 2 ) = (dt) 2 -a(t) 2 (dx 1 ) 2 . (1) 

The formulation of the Newtonian problem in this metrics proposes a choice of physical variables 
and coordinates. The modern cosmology uses two choices of such variables: the conformal time (r]) 



and coordinate distance (x 1 ) with the interval 

(ds 2 ) = a( V ) 2 Udr]) 2 - (dx 1 ) 2 



(2) 



and the Friedmann time (t) and distance X % 



ax 



in terms of which the interval (0) takes the form 



(ds 2 ) = (dt) 2 - dX l - H(t)X l dt 



(3) 



where H(t) = d(t)/a(t) is the Hubble parameter, we have used here the formula of differential calculus 
adx = d(ax) — xda. Both the sets are mathematically equivalent 1 . 

In terms of the Friedmann variables X % = ax 1 the Newton action in the space with the interval 
(J3j) takes the form 



to 

Sa = J dt 
ti 



Pi(X l - HX l ) - 



P 2 



2m 



(4) 



1 Note that a(t) = a(r](t)) is connected with a(rj) by equation dt — a(rj)dr]. 



The equations of motion 

X* - HX l = — . P l + HP 1 = 

have the simplest solution X 1 where constant. 



m ' 



(5) 



The transition to the conformal variables gives the action 2 



S A = J drj 

VI 



Pi 



dx l 



Pi 



drj 2m a(r]) 



(6) 



for a particle with the running mass moa(r)) (with the present-day value a(rj ) — Oo = 1). This 
transition is just the main idea of our paper to separate the Hubble velocity field 



H,. 



1 dR 
~R~dt 



R=^xt + xl + xi) 



(7) 



from possible Newtonian motion. Substituting R = ar, where r = Jx\ + x\ + x|, in the definition 
of the total Hubble flow (J2J) we get this Hubble flow in the following form 



Ida 1 dr 
Htot = -— + -— = H + AH. 

a dt r dt 



(8) 



One can see that the total Hubble flow differs from the classical one H = a/a by the value of the 
local flow 



AH 



1 dr 
r dt 



(9) 



III. NEWTONIAN MOTION IN AN EXPANDING UNIVERSE 



Let us consider the action 



t 



Sa = J dt 

ti 



P 2 a 

Pi(X l -HX l ) *- + - 

2mo R 



(10) 



where a = Mom^G is a constant of a Newtonian interaction of a galaxy with a mass mo in a 
gravitational field of a central mass Mq. Action (fTU|) for radial momentum Pr and orbital moment 
Pq in the cylindrical coordinates 



X l = Rcos6, X 2 = RsinO, X 3 = 



2 The action (@J and JHJ differ from the one (7.4) in monograph 0, that is not compatible with quantum field theory 
on the conformal flat metric (2) 



takes the form 

?2 



S A = fdt 



;i2) 



ti 

To separate the Newtonian motion from the Hubble velocity field, we use to the conformal variables 
Pr = -Pro(^), t = R/a(t), drj = dt/a(t). In terms of these variables the action takes the form 



S A = J dr] 



V ° r r? P 2 a 

p y + p e e' Pi tl + _ 

2moa(r]) 2moa(rj)r 2 r 



VI 



(13) 



IV. THE CAPTURE OF GALAXIES BY THE CENTRAL GRAVITATIONAL FIELD 

The energy of a particle with the running mass m(r]) = a(t])mQ described by the action (fT^j) 

E{rj) = ^ + — — (14) 

2moa(r]) 2moa(r])r 2 r 

is not conserved in the contrast to the energy of particle with a constant mass in the Newtonian 

mechanics j5|. In our case (|14j) . if the scale factor a(r]) increases, the energy (|14|) runs from its 

positive values to negatives ones. There is a moment of a time 7] = r]i when the energy ()14|) is equal 

to zero: 

J^L s W^A_^ = 0i (15) 

m a(?7/) 2 

where r\ = p r {vi)/ m i r i is radial initial velocity, mj = m a(rij), rj = r(r]i) are the initial conformal 
mass and coordinate distance, and 



vi = — — , wi = W—^— (16) 
rrijrj y mjrj 

are the orbital velocity and Newtonian one, respectively. It is known that the change of a sign of 
the energy means the change of an unrestricted motion of a particle by a finite motion in the central 
field. Therefore, the time rji can be treated as the time of the capture of a particle (cosmic object) 
by the central gravitational field. 

If the initial radial velocity is too equal to zero r\ = 0, the zero-energy constraint ()15j) 

vj = 2wj (17) 

becomes the equation for the initial data moa(r]j)r(r]j) = mjrj. The solution of this equation 
rrijrj = Pq j i2a) can give a orbital velocity for all trajectories of the captured cosmic objects 

2a 

Vj = — = constant (18) 
Pff 



The fact of the universality of the orbital velocity (jl8[) for all ellipsoidal trajectories (due to the zero 
energy initial data of the formation of a local universe) gives us a possibilityof expanding both the 
o_ observational data on the law of the constant orbital velocity BflBQ and the anizotropy 
of the local velocity field 3, 4 . 

nn 

The anisotropy of the local velocity field AH = [jj, |2( is not compatible with the class of the 
isotropic circular trajectories r(no) = ro,r' = r" = with the equation of motion 



vl = wl lv = w = — ) (19) 

V m r V m o r o J 



known as the "virial theorem" [3, 18, 9], where the initial radius rj can be identified with the obser- 
vational radius 3 . If trajectories belong to a class of ellipsoidal ones, the initial radius rj does not 
coincide with the observational radius. 



V. THE EXACT SOLUTION 



Let us consider a solution of the Kepler problem, i.e., a motion of an object in the central field in 
the expanding Universe for the rigid state when the densities of energy and pressure are equal. The 
cosmic scale factor a can be written as Q| 



a(r]f 



which describes the recent Supernova data 



[1 + 2H ( V - no)] 

□ 



11 



12 



in the relative units 



0. 



(20) 

Equation (j2*Uj) can 



be considered as a change of the evolution parameter n —* a {v)- m terms of the new variables 

y = r/r , p y = p r /m , v = P 9 /r m , 
where r(no) — r o is the present-day data, action (jlHj) takes the form 



S A = r m J da j p y ^ + v 



de 

' da 



1 

c 



V P 2 y + vl/y 2 aw, 



(21) 



where wo = ^Ja/m r Q , Co = H ro are the Newtonian velocity and the Hubble one, respectively, and 
aj = 1/(1 + zj) is determined with the redshift Zi at the moment of formation 77 = ni- The total 
energy of the system (|T4*|) in terms of the new variable takes the form 



E(a) 
r m Q 



1 

Co 



P 2 y + vl/v 2 _ awl 
2 y 



(22) 



3 Remind that the hypothesis of the Cold Dark Matter was proposed in |7], |8| to explain the law of the constant 
orbital velocity in the class of the circular trajectories in the nonexpanding Universe where l|19|l is valid. 



It is easy to see that vq is a constant of the motion dvo/da = 0. The equations of the radial motion 



dy 

can be written in the Lagrangian form 



c 



dp y Vq awl 



da y 3 y 2 



(23) 



m 2 




)'- 


a 






> y 3 


y 2 _ 



da 2 

At the present-day time a = ao = 1, yo — 1 this equation determines the second derivative 



(24) 



<fy 

da 2 



a=l 



Wq 
C 



Wq 



- 1 



(25) 



in terms of three velocities u>o,^ojCo, and it allows us to choose any relations between wq and vq, in 
particular v 2 = 2w 2 , in contrast to the circular trajectory where Vq = w%. 

The general solution of this equation (j2I| can be obtained in the parametrical form: 

N 2 (r) 



a t 



cr 



y( T ) = C2 r 2 / 3 iV(r) 



r 2/3 ]V ( r ) ' 

where r is the parameter of solution with the initial date r = r(i] ), and 



(26) 



N(t) = aiU 2 {T)+/3 1 U{T)V(T)+y 1 V 2 (T), 
N 2 (r) 



r dN(ri + 2 i ±4 r 2 iV 2 (r)+^A, 
dr 3 I 



A = 4a l7l - /? 2 : 

f v \ f 3c 
ci = — 



1/3 



wo/ Uw / 2a;A 1 /2' 
w \ {4w \ 1/3 1 



vWq/ \ 3co / a; A 1 / 2 

Three constants a.\, 71 = const can be found from the following system of three equations: 

1 



(27) 
(28) 
(29) 
(30) 

(31) 



lr=0 



a\ = aj = . 

1*7=77/ 1 l + Zl 



I ri 1 1 ^ 

y r=0 = — = 1 + z h T 

r da 



0: 



(32) 



r=0 



here for the upper sign (restricted solution at infinity, r = +00) in ()28j) 

C/(r) = J 1/3 (r), V(t) = y 1/3 (r), and u = - 



(33) 



where Ji/s(t) and Yi^r) are the Bessel functions of the first and second (or Niemann function) 
kind, while for the lower sign (unrestricted solution at infinity, r = +00) in (J28j) 



U{t) = J 1/3 (r), V{r) = K 1/3 {t), and 



(34) 



where I\/%{t) and Km^t) are the modified Bessel functions of the first and second (or MacDonald 



function) kind (see, e.g. |l4j). 

A solution of the equation of motion following from the action (fTHJ) is given in Fig. ^and Fig. 121 
We use this solution for construction of two plots. 

Figure 01 gives the values of the correction © to the Friedmann Hubble flow resulting from taking 
into account eq. ©• It is clearly seen that the corrections are dumped with time and have a 
quasi-periodical character. 

In Fig. 13] the angular distribution of the Hubble flow correction, as given by @, is presented. 
The correction is anisotropic. We consider the 2-dimensional case while Karachentsev's anisotropy is 
observed in 3-dimensions. Nevertheless, our 2-dimensional analysis allows one to see the anisotropy 
of the Hubble flow and to estimate the order of the magnitude of the anisotropy. 

We have considered the case when the formation of a galaxy began from the zero-energy state 
(the initial data E = 0, and velocity y f = 0). These data correspond to the relation = 2w^. 

VI. CONCLUSION 

Our paper was motivated by the finding that in the local Universe the velocity field is anisotropic P, 
This effect is difficult for explanation. The only possible suggestion, but rejected by Karachentsev, 
was rotation [4] . We are trying to find the origin of this anisotropy. In ordered to do this, we consider 
the general uniform expansion of the Universe. Since this is the nearby (less than 8 Mpc) part of 
the Universe around us, we use the Newtonian approach. We studied the motion of the test massive 
particle in the central gravitational field on the background of cosmic evolution of the type of FLRW 
space-time with uniform expansion. We assume the rigid state of the matter when densities of 
energy and pressure are equal and which corresponds to conformal cosmology compatibles with 



Supernova data . We obtained the exact solution of the above-mentioned Kepler problem, to find 
the difference between the uniform Hubble flow and our case. We have shown that this difference was 
anisotropic. In such a way we explained the anisotropy of the local velocity field by the Newtonian 
motion of galaxies in the central field. Of course, our 2-D consideration shows a possible mechanism 
of the observed 3-D anisotropy. 

Having the solution of the Kepler problem we admit the rotation of galaxies around the center of 
the local Universe. This local Universe must be regarded as the Local Group of galaxies. In such 
a way, we support the picture in which galaxies rotate around the center of the Local Group in the 



class of ellipsoidal trajectories. 
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Figure 1: Graph of function y(r) = a{r)r{r) /tq (|26() 
at vl/cl = 0.2 and wl/cl = 0.1. 
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Figure 2: Graph of functions product a(r)y(r) (|26|) 
at ug/c^ = 0.2 and wg/c§ = 0.1 for -3 < a(r) < 14. 




Figure 3: Graph of function AH in units Hq at Figure 4: Graph of function AH in units Hq at 
ug/c§ = 0.2 and wl/c 2 = 0.1. wg/c^ = 0.2 and wl/c 2 = 0.1 for 0.8 < a(r) < 14. 



